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ABSTRACT 



Aims. We investigate the instability of toroidal magnetic fields resulting from the action of z-dependent differential rotation on a 
given axial field B° in a cylindrical enclosure where, in particular, the helicity of the resulting nonaxisymmetric flow is of interest. We 
probe the idea that helicity is related to the external field and the differential rotation asH on B° B° Qij. 

Methods. We conduct isothermal magnetohydrodynamic simulations of a quasi-incompressible medium with finite viscosity and 
conductivity in a perfectly conducting container, and analyze both the kinematic and current helicity of the resulting field by regarding 
the nonaxisymmetric parts of the field as fluctuations. 

Results. The observed instability leads to a nonaxisymmetric solution with dominating mode m = I. With the onset of instability, 
both kinematic and current helicity are produced which fulfill the suggested relation H oc B'l B'j Qj j. Obviously, differential rotation 
dQ/dz only needs an axial field S? to produce significant helicity. Any regular time-dependency of the helicity could not be found. 
The resulting axial a-effect a- is mainly due to the current helicity, the characteristic time scale between both the values is of the 
order of the rotation time. If the axial field is switched off, then the helicity and the a-effect disappear, and a dynamo is not observed. 
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1. Introduction 

Helicity in rotating turbulence plays a basic role in dynamo 
theory for the generatio n of large-scale cosmic magnetic field s 
jKrause & Radler 19801: iBrandenburg & Subramanian 20051) . 
Usually stratification of density or turbulence itself is needed 
to generate helicity. We demonstrate that, even without strat- 
ification, helicity can exist due to the common action of 
magnetic fields and differential rotation. Whereas differential 
rotation depending on latitude is well-known from observations 
of star spots, differential rotation (of the interior of stars) 
depending on the z-coordinate only, appears more rarely. It 
is predicted and observable, for instance, in the tachocline of 
the sun (Rudiger & Kitchatinov 1997; Kitchatinov & Rudiger 
2005), and influences th e internal rotation of massive stars 
jMaeder & Mevnet 200 3\ It may open the possibility of study- 
ing nonaxisymmetric (large-scale) structures by shear-driven 
magnetohydrodynamic (MHD) instability o f laminar flows 
jBraifliwaite & Spruit 20041: iBraifliwaite 2006h . 
Differential rotation transforms poloidal field components into 
toroidal components, which due to the Tayler instability (TI), 
become unstable if a critical amplitude is exceeded (Vandakurov 
1972). Tayler (1973) showed that a magnetic field becomes 
unstable against nonaxisymmetric perturbations if the condition 



smoothing action of the differential rotation. In a mean-field ap- 
proach, due to the relation 



(2) 



(1) 



helicity 'H - TYkin - 'Hcun often indicates an existing a-effect. 
•Kkin and 'Hcun are the kinematic and current helicity, 

-Kkin = <h' ■ l-OtH'), •T/eurr = —{B' ■ lOtB'), (3) 

/'OP 

respectively and T^on- is the correlation time of the field pattern. 
Nonvanishing helicity could indeed be an indication of the ex- 
istence of a resulting a-effect on the basis of the TI. In case of 
the presence of an external field B'^ and with a large magnetic 
Reynolds number, a dominating contribution from magnetic- 
field fluctuations to the a-effect can arise (Pouquet, Frisch & 
Leorat 1976: Brandenburg & Subramanian 2007). 
Keeping in mind the pseudoscalar nature of helicity, a possible 
relation between helicity and both the external field and differ- 
ential rotation would be 

9{ ^ B'l B^'j Qij, (4) 

where Z?'' means the external field components and Qjj the gra- 
dient tensor of the basic rotation. The sign of 'K in Eq. (|4]i does 
not depend on the sign of the magnetic field, but it depends on 
the sign of the shear Qij. Consequently, the dynamo number 



where cylindrical coordinates (R,(p,z) are used, is violated. 
Constant fields or fields B^ oc like those produced by z- 
dependent differential rotation from an axial field, are candidates 
for this instability. The most unstable Fourier mode for such a 
configuration is the nonaxisymmetric mode m - 1 despite the 
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would now have a positive-definite sign. Here L is the character- 
istic scale of the cylinder domain and rjj ^ Tcoit(m'^)- The sign 
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of the dynamo number determines important properties of pos- 
sible dynamo processes. E.g., a simple disk dynamo oscillates 
for a positive dynamo number and it is stationary for a nega- 
tive dynamo number Also, the magnetic Reynolds number of 
the (T-eflFect itself is important, that is, with Eq. Q including the 
magnetic contribution. 



(77T does not include magnetic contributions, Vainshtein & 
Kitchatinov 1983). We shall see that after the onset of TI, this 
quantity is of the order of unity when the averages are taken as 
an integral over the azimuthal direction. 

In the following, the details of the model used are explained. 
Section|3]describes the onset of the TI in a cylindrical enclosure 
and the resulting field structure. In Sect. H] it is shown that - if 
the axial field is strong enough - the numerical simulations are 
in good agreement with relation|4]in the case of z-dependent dif- 
ferential rotation, where it simplifies to 'K oc B^fdQIAz and only 
a (strong enough) z-component of the external field is needed. 
In Sect. |5]we show that the generated helicity is connected with 
an a-effect, which is also proportional to the gradient of the dif- 
ferential rotation. 



To keep the system nearly incompressible with a code for com- 
pressible media, the value of Q{z) is chosen in a way that the 
resulting maximum of the (meridional) velocity never exceeds 
16% of the speed of sound, which was set to Cs = 8. The ini- 
tial density is set to po - 1, fluctuations are of the order of 
10 The spatial resolution is = 96x96x48. Test runs with 
= 128x128x96 deliver only slight differences with smoother 
fields. 

The initial magnetic field consists of the external time- 
independent homogeneous field = (0, 0, B?) applied only in- 
side the cylinder It is twisted into a strong toroidal field by the 
differential rotation. The strength of the generated field depends 
on B*', varied between 5° = 0.01 and 5° = 0.1, and the gradient 
of the differential rotation AQjAz, where the latter was fixed to 
|di3/dz| = 1 for all presented calculations. Based on the radius 
(set to L = 1) and the velocity si z - \ {U - QL - 1), and 
the viscosity (v = 0.03), the Reynolds number has a value of 
Re = 33. In the following text, system rotation always refers to 
the rotation on top of the cylinder and time is given in units of ro- 
tation time. The magnetic Prandtl number Pm is varied between 
Pm =10 and Pm - 30 with no qualitatively change of the in- 
stability. In the following analysis, we concentrate on Pm =15. 
We choose a large Prandtl number to encourage dynamo action, 
although we observed no dynamo. 



2. Model 

We consider a differentially-rotating cylinder with radius L and 
height L embedded in a box with cartesian grid of side length 
2.5L and height L (see Fig. [T]i. The cylinder radius is 20% less 
than the box side length, which appeared to be a good compro- 
mise between wasted computing power in the corner regions and 
decreasing influence from the geom etry of the box. We calcu- 
lated using the PENCIL code (Bran denburg & Dobler 20 02). a 
high-order finite difference code (sixth order in space and third 
order in time). Inside the cylinder domain, we solve both mo- 
mentum and induction equations. Outside the cylinder, we kept 
the initial velocity to sustain the differential rotation via this no- 
slip conditions on the rim of the cylinder. On the top and the bot- 
tom, we apply stress-free conditions for the flow. For the mag- 
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Fig. 1. The considered domain is a cylinder inside a box with 
Cartesian grid of side length 2.5L and height L. The full set of 
equations is solved only in the cylindrical domain. 



netic field, perfect conductor boundary conditions are applied on 
all box boundaries. Additionally, in the region outside the cylin- 
der, magnetic diffusivity r\ is enhanced to a value ten times larger 
than inside the cylinder 



3. Instability 

If becomes strong enough, the TI occurs and leads to a grow- 
ing nonaxisymmetric field. The largest nonaxisymmetric mode 
is m - 1. By "strong enough" we mean that not only does 
the magnitude of reach high enough values, but also a cer- 
tain threshold B^jB^. (with B^jB^} of the order of Pm) needs 
to be crossed, which is different from the case of pure toroidal 
fields (Tayler 1973; see Riidiger et al. 2007). This means that 
an additional poloidal field component suppresses the instabil- 
ity. For Pm = 15, the instability sets in at a Hartmann number 
of Ha = B'Rj yfJJoQvrj - 130, where B* » 1 means the maximal 
value of the generated toroidal field B^. To produce such a strong 
toroidal component an external field B° > 0.04 is needed. The 
instability did not occur forB'^ > 0.08. 

The nonaxisymmetric structures appear first near the axis in the 
lower part of the cylinder where velocity magnitude is small, 
and grows to the steady state (Fig.|2]i. For Br and B , the m - 1 
mode becomes the largest one, in B^ the axisymmetric mode 
remains dominant due to the permanently reproduced axisym- 
metric field. In the nonlinear regime, before the steady state is 
reached, higher modes also appear Nonetheless, their magni- 
tudes stay below 10% of the m - 1 mode. One exception is the 
m = 4 mode, caused by the box geometry. During the initial 
phase, this mode is already present and influences the magnetic 
field outside the cylinder near the box boundaries. However, the 
field inside the cylinder is much less affected, and the m = 4 
mode does not seem to influence the nature and onset of the in- 
stability at all. Additionally, during the growth phase of all other 
nonaxisymmetric modes, the m = 4 remains nearly unchanged 
and does not play an extraordinary role in the final state. Modes 
higher than m = 4 are not plotted in Fig. |2] 
The pattern of the m - \ mode (Fig. |3]l possesses an az- 
imuthal drift velocity relative to the system rotation. It rotates 
with 4.5 rotation periods of the cylinder, independent of the 
strength of the externally applied field. The TI in a Taylor- 
Couett e system exhibits th e same characteristics of the unstable 
mode jRiidiger et al. 2007h . 
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Fig. 2. The power spectrum of the Fourier modes of for Pm - 
15 and = 0.05. Onset of the instability near f = 28 where the 
dominating nonaxisymmetric mode m - \ starts to grow. Mode 
w = 4 is nearly unaffected by the instability. 
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Fig. 3. The m - 1 pattern in the magnetic field for Pm =15 
and - 0.05. The pattern drifts relative to the systems rotation 
with 4.5 times the rotation period of the cylinder. The cylinder 
boundary is marked by the dashed black-white line. 



4. Generation of helicity 

To test the relation described by Eq. (|4]l in a simple setup, we re- 
stricted differential rotation to depend only on the z-direction. 
For this case, the helicity 'H should depend on the sign of dQ/dz, 
but not on the sign of the external field and scale with the squared 
value of the latter, i.e. "H oc B^}^ dQ/dz. During the onset of the 
instability, the deviation of the flow structure from the original 
toroidal field is comparable to that of the magnetic field. Also, 
here the mode m = 1 is the largest nonaxisymmetric mode and 



' The resulting non-conservative centrifugal force drives a merid- 
ional flow that, however, remains small (see Fig.O. 
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Fig. 4. The external field dependence on the energy of velocity 
and magnetic field fluctuations at /? = 0.5 averaged in vertical 
direction for 0.35 < z < 0.65. 

the higher modes appear with less energy. 
For the definition of fluctuations, we average velocity and mag- 
netic field along the azimuthal direction, i.e. 



u d<p. 



Bd4>, 



(7) 



and regard the nonaxisymmetric parts as fluctuations of the 
fields. For constant dQ/dz, and with a homogeneous external 
field, the fluctuations should be constant in the vertical direc- 
tion. Due to the finite height of the cylinder, boundary effects 
near the top and bottom induce additional flow disturbances and 
changes in the structure of the magnetic field in these regions. 
The independence is roughly preserved in the central region be- 
tween z - 0.35 and z - 0.65. In Fig. |4]we show the energies 
of the fluctuations of velocity and magnetic field averaged over 
this central region, and their dependence on the external field B?. 
The values are taken in radial direction ?AR - 0.5. The magnetic 
energy always exceeds the kinetic energy. The ratio varies be- 
tween 15 and 30. Thus, the instability is dominated by the mag- 
netic field. This is also reflected in the helicities. We again use 
azimuthal-averaged quantities. 



n 



kin ■ 



rot u' dip, 'Ha 



1 



B' -mlB'dip. (8) 



Both are zero as long as no instability exists, with proceeding in- 
stability and growing nonaxisymmetric modes the helicity quan- 
tities also become unequal to zero. An example plot of "Hkin for 
R - 0.5 and z - 0.5 is shown in Fig. |5] After a transition, a 
steady-state value is reached for both "TYkin and "TYcun . The cur- 
rent helicity is approximately 40 times larger than the kinetic 
helicity in the example shown in the center of the cylinder. In 
the vertical direction we observed the same behavior as for the 
fluctuations of flow and magnetic field. Near the top and bot- 
tom, the helicities are stronger than those in the region between 
z = 0.35 and z = 0.65. For a positive gradient in angular velocity 
TYkin is negative and 'Kcun positive. With the changing sign of 
dQ/dz, the sign of "Hkin and "Hcun also changes, whereas the sign 
of B? has no influence. And, as expected, the absolute value of 
'^kin is proportional to the squared value of the external field. For 
Pm = 15 the relation between B? and TYkin is shown in Fig.|6] 
The fitted parabola does not exactly cross the origin. This shift 
is attributed to the instability requiring a certain value of B? (i.e. 
B[? = 0.04 for Pm = 15) to set in with a certain value for "Hkin 
larger than zero. 
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Fig. 5. Kinematic helicity 'Kkin and current helicity 'H^^^ 2lVR - 
0.5 and z - 0.5 for Pm = 15: both have different signs whereas 
signs of current helicity and gradient of angular velocity df2/dz 
are the same. The external field values are B^l - 0.05 for negative 
and B° - 0.06 for positive gradient. Note the different size of 
-Hkin and 'Hcun- 
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Fig. 6. Kinematic helicity "Hidn as function of the external field 
with fitted parabola. Values of TYidn are taken at R - 0.5 and 
averaged along vertical direction for 0.35 < z < 0.65. 



5. Alpha effect 

Equation (|2]i suggests that the changing sign of the helicities by 
changing the direction of dQ/dz will also appear for the term 
a,, in the a-tensor responsible for the regeneration of the mag- 
netic field. Indeed, calculating a., = {u' x B'),/B^^ using the z- 
component of the electromotive force, gives values for a,, with 
the same sign as dQ/dz for z < 0.65 (see Fig.|7J. Near the top 
boundary, the disturbances already discussed for the helicities 
change the values of dramatically to the opposite sign. We 
excluded this region and took mean values of a^r in the same 
range 0.35 < z < 0.65 used for the helicities and find coeffi- 
cients between a^, = 0.038 for B° = 0.04 and a,, = 0.067 for 
B? - 0.08. The correlation time derived from expression dU us- 
ing both helicities, Tcorr - 3a,.j/('?Ycurr - 'Hkin), is of the order of 
the system rotation time. Note that the a-effect here is mainly 
due to the current helicity and the product a^/Hcun is positive. If 
the external field is switched off, all magnetic field modes decay 
in the investigated parameter region of low Reynolds numbers 
and magnetic Prandtl numbers of the order of ten. The helicities 
drop down to zero nearly immediately. Despite the high mag- 
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Fig. 7. Coefficient of the a-effect a^z aXR - 0.5 for positive and 
negative dQ/dz. Sign of a~, is the same as 'Hcun and dQ/dz. 



netic Prandtl number, we do not observe any dyna mo action, 
as found in comparable geometry (Br aithwaite 200a) . Also, the 
reported cyclic behavior of the instability did not occur in our 
simulations. 



6. Conclusions 

We have demonstrated how the instability of toroidal magnetic 
fields leads to helicity generation in the nonaxisymmetric parts 
of a flow without density stratification. The suggested relation 
between helicity, external field, and differential rotation, 'H oc 
B^l B^j Qjj, is rather well fulfilled in the case of z-dependent dif- 
ferential rotation. Both the kinematic helicity and current helic- 
ity depend on the squared value of the z-component of the exter- 
nal magnetic field Z?*?, and scale linearly with the gradient of the 
differential rotation. Also, the magnetically dominated a-effect 
depends on the direction of dQ/dz, that is a,- holds the same 
sign as dQ/dz. 

The realized model is too simple to estimate consequences of an 
a-effect based on this new kind of helicity production for envi- 
ronments without density stratification and possibly for new dy- 
namo models. In a next step, therefore, we would like to check 
the relation (HI when differential rotation depending on the dis- 
tance from the rotation axis is present as well. With an appro- 
priate flow it becomes 'Kkin ^ B^j^B'^dQ/dR. In this case, helicity 
generation should be observable only if both components of the 
external field are unequal to zero. The problem with this constel- 
lation is the changing sign of the product of both field compo- 
nents in one hemisphere in the simplest configuration, a dipole- 
like field. 



References 

Braithwaite, J. 2006, A&A, 449, 451 

Braithwaite, J., & Spruit, H.C. 2004, Nature, 43, 819 

Brandenburg, A., & Dobler, W. 2002, Comp. Phys. Comm., 147, 471 

Brandenburg, A., & Subramanian, K. 2005, Physics Reports, 417(1-4), 1 

Brandenburg, A., & Subramanian, K. 2007, A&A, in prep. 

Kitchatinov, L.L., & Riidiger, G. 2005, Astron. Nachr, 326, 379 

Ki'ause,F. , & Radler, K.-H., Mean-field magnetohydrodynamics and dynamo 

theory (Akademie-Verlag, Berlin, 1980) 

Maeder, A., & Meynet, G. 2003, A&A, 411, 543 

Pouquet, A., Frisch, U., & Leorat, J. 1976, JFM, 77, 321 

Riidiger, G., & Kitchatinov, L.L. 1997, Astron. Nachr, 318, 273 

Riidiger G., Hollerbach, R., Schultz, M., & Elstner D. 2007, MNRAS, acc. 

Tayler, R.J. 1973, MNRAS, 161, 365 

Vainshtein, S.I., & Kitchatinov, L.L. 1983, GAFD, 24, 273 

Vandakurov, Yu.V. 1972, SvA, 16, 265 



